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containing inverse powers of the Ricci scalar varied with respect to the metric. While 
we can fit the supernova data well, we confirm the a oc t 1 / 2 behaviour at medium to 
high redshifts reported elsewhere and argue that the easiest way to show that this 
class of models are inconsistent with the data is by considering the thickness of the 
last scattering surface. For the best fit parameters to the supernova data, the simplest 
1/R model gives rise to a last scattering surface of thickness Az ~ 530, inconsistent 
with observations. 
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1. Introduction 

Observations of type la supernovae still consistently suggest that the universe is 
accelerating [H [21 [3] (see [I] for a detailed list of references). This conclusion seems 
to be supported by the fact that observations of the CMB tell us that the universe is 
spatially flat [5] and that the matter which seems to be responsible for galaxy clustering 
is not enough to account for this flatness 0,[7]. Furthermore, the measurement of space- 
time geometry achieved by the detection of the imprint of the waves in the primordial 
plasma in the galaxy correlation function [8] has caused problems for models where the 
supernova data are explained by the dimming of photons P, QUI tH]- 

There have been a number of different approaches trying to explain the mystery 
of the acceleration of the universe. One such approach consists of a class of theories 
obtained when the Einstein-Hilbert Lagrangian is modified by hand and the Ricci scalar 
R is replaced with some function f(R) [121 [T3]. These theories are phenomeno logical 
in as much as it is not clear what the underlying theory that gives rise to them would 
be. The motivation is that inverse powers of the curvature in the Lagrangian will 
give vacuum solutions which are not Minkowski but which rather lead to late time 
acceleration. The simplest theory with inverse powers of the curvature in the Lagrangian 
is 

S=^ff d'x^-g (R-Ji)+J d 4 x^- 9 C M (1) 
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where \i is a mass scale that must be fitted to the data. Variation of this action with 
respect to the metric gives new field equations 

+ jp) - 2 y- - jp) Rg ^ + ^v a v a - v M v,];r 2 = ^ (2) 

from which it can be seen that for cosmological solutions there will exist a vacuum 
solution with H ~ \i. Such f(R) models can therefore give rise to late time acceleration 
which could be responsible for the apparent dark energy. 

The simplest models of this nature seem to be at odds with solar system tests 
of gravity, at best containing a light degree of freedom [H] and at worst possessing 
instabilities [H2, [16] but one might imagine that the Lagrangian (TjQ) is some effective 
limit valid on very large scales and that some new physics on short distances could 
change the theory. There are also indications that some f(R) theories may be safe in 
some regions of parameter space [T7J [HI EH] through a process rather similar to the 
chameleon mechanism [20J although that does not appear to be the case for the ones 
studied here. 

It should be noted that variation of the normal Einstein-Hilbert Lagrangian with 
respect to the metric or alternatively the Christoffel symbols leads to the same field 
equations for gravity whereas in this class of modified Lagrangians the same is not 
true. Field equations obtained using the latter Palatini approach will yield different 
cosmologies and solar system constraints [211 1221 [23]. In this work we will restrict 
ourselves to equations of motion obtained by varying the Lagrangian with respect to 
the metric. 

The solution of the field equations for a cosmological background lead to 
cosmological equations with higher derivative terms, for example for the field equations 
(J2|) the tt Friedman equation for a spatially flat universe becomes 

3 #2 / 2HIJ + lhH 2jj ^2 q H A = PM (3) 

12(H + 2H 2 ) 3 V ) M 2 Pl y ' 

which means that there are more degrees of freedom in the space of solutions than for 
Einstein gravity with a cosmological constant or the standard DGP model. This space 
of solutions needs to be compared with the data. 

If type la supernovae are good examples of standard candles as is thought, they 
can in principle trace out the Hubble diagram in an unambiguous fashion. There is, 
however, a need for some caution as there are presently ambiguities in the way that 
one can analyse the data [24, [25] and also in which supernovae should be included in 
samples to be used for cosmology [H [2], [3] . In certain situations, these ambiguities can 
lead to different predictions with regards to which dark energy models are favoured over 
others [251 [23 S] . 

Despite this, all of the supernova surveys seem to agree on some basic facts, namely 
that the universe is accelerating and that ACDM, i.e. a universe composed of matter, 
radiation and a constant energy density which does not change over time, fits the data 
rather well, and better than many alternatives motivated by specific physical models. 
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Expanding the surveyed redshift range is of great importance to investigate the 
nature of dark energy and to find out in particular whether it is a cosmological constant 
or not. Thus, the recent compilation of supernova data in [2], including 23 SNIa at 
z > 1 from an HST/ACS is a very interesting data set for model comparisons. Riess 
et al have included data from various sources and re-fitted the light-curves with the 
MLCS2k2 technique. 

We have used the gold set of 182 SNe la from Riess et al [2], taking into account 
the additional redshift error discussed there in the case of the high redshift supernovae 
for which the redshifts were determined from broad features in the spectra. 

Furthermore we have added to this data set the position of the acoustic peak in the 
SDSS galaxy survey [8] and also the CMB shift parameter [28]. We direct the reader to 
[4] for more details and discussion on measuring the expansion history with these data 
sets. 

This data leads to a set of angular distances and luminosity distances which need 
to be explained by any successful dark energy model. In this work we aim to find out if 
it is possible to fit this data using the theories of modified gravity introduced above. We 
also need to verify that the best fit models do not at the same time change conditions in 
the early universe such that the CMB and baryon oscillation data is no longer consistent 
with the physics of that epoch. 

In other words, in order for the models to be acceptable, there has to be a circular 
self consistency between the theory and the data. It has been observed that if such 1/R 
models give rise to acceleration at late times, then they also lead to expansion at early 
times which mimics, in terms of effective equation of state, the expansion during the 
radiation dominated era [30l l3T] (independent of whether or not there actually is any 
radiation in the universe). One question which needs to be answered therefore is what 
is the density of this effective radiation component for model parameters which also fit 
the observed luminosity and angular distances in the data. Having obtained the value 
of that effective density, we need to see if it is consistent with the observed properties 
of the CMB. If not, then the theory may appear to be able to fit the existing data but 
the data is no longer consistent with the early universe physics implied by the theory, 
and the model is ruled out. Finding out if it is possible to evade this situation, and fit 
all the data simultaneously is the subject of this paper. 

2. Obtaining solutions for f(R) theories. 

There are a large number of different f(R) theories with inverse powers of R that one 
might consider but we will look at the two simplest models that can be written down, 
namely f(R) = R- tf/R and f(R)=R- fi 6 /R 2 . 

In order to test the f(R) theories, we need H(z), the Hubble expansion rate as 
a function of redshift. This is rather difficult to obtain from higher order Friedman 
equations such as (j3J). We choose to perform a conformal transformation which allows 
us to treat the problem as one of a scalar field a with a potential V(o~) evolving in an 
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FRW universe where the matter is endowed with a non-standard gravitational coupling 
set by a. 

Having found the solutions in the Einstein frame for the scale factor, Hubble 
parameter and the scalar field, we will transform back to the matter frame where 
we will integrate H(z) and compare it with the data. 

For an action of the form f(R) with Of /OR > 0, we can perform a conformal 
transformation to the Einstein frame [29|. The conformal transformation is written 

9^=P9^, ^=P = ^[VV^^ S ) ( 4 ) 

we will continue to use both p and a even in the same equations to keep expressions 
compact despite the fact that they can be used interchangeably. The equations can now 
be written in the more familiar form 

R, u - l -R~g, v = (V^)V„o- - l -g^{V a a)V p a - V{a)g^ + (5) 

where the potential and the energy-momentum tensor in the Einstein frame are given 
by [29] 

V = {si9n) (r^-A (6) 
2\df/dR\ V dR J J { ' 

f, v = ^ (7) 
V 

with (sign) = J^/ | The potential V has mass dimension 2 because we have chosen 
to work with a dimensionless scalar. We denote quantities in the Einstein frame with a 
tilde, for example the time coordinate in the Einstein frame Robert son- Walker metric 
is dt = y/pdt and the scale factor a(t) = v /pa(t). The equations which need to be solved 
simultaneously are the equation of motion for a and the Friedman equation. The latter 
comes from the time-time component of the transformed Einstein equations 

a* -l,- + (8, 

where a prime denotes differentiation with respect to t, the density pu = Pm/p 2 and 
H = a! /a. The equation of motion can be obtained from the time component of the 
divergence or from the covariant derivative of the stress energy tensor T 

Equations (jSJ) and (jSJ) are the ones that we need to solve to get the evolution of the 
universe in these models. In order to relate the Hubble parameter H in the matter 
frame and H in the conformal frame, one must use 

do) 

| the frame where the matter is not coupled to gravity via a, sometimes called the Jordan frame or 
more recently the string frame 
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To compare these models with the data, we need H(z)/H in the matter frame, or rather 
J H (gy Ho ■ There are more free parameters than in ACDM so we restrict ourselves to the 
case k = 0. For a flat universe, our free parameters are VIm as defined in the original 
matter frame, the mass scale /i and the values today of a and a'. To write these in 
dimensionless form, let us introduce the time-variable r = Hot and the parameters 



#o R , I da 

fj, V6 0T 



where the subscript denotes the values today at redshift z = 0. If we also introduce 
the dimensionless potential U(p) so that 

V(a) = /2 2 U(a) (12) 

then the rr Friedman equation becomes 

| (13) 

Hi 6 \ot J 3 a 2 y/pap \ a J 

while the equation of the motion for the scalar field is written 

av + 3 ^ + ^» _/3t>^c*y = . (14) 



dr 2 H dr a 2 da \ 2 y/pop \ a 

By solving ffTBI) and f[T4"|) numerically, we can recover the evolution of the Hubble 
expansion in the Einstein frame and can convert this back to the matter frame using 
o~(t). One of the parameters can be written in terms of the others using (TlBl) with z = 



2_ PlPUjPo) n ,x 

3(2p -p (3-n M (3) { ) 

We therefore choose to label the parameter space of solutions of H(z)/H with the three 
variables Qm, o"o and da/dr\ . 



3. The data vs. f(R) =R- fi 4 /R and f(R) = R- ^ / R 2 

In this section we will see if it is possible to fit the existing data under the assumption 
that the effective dark energy density is not significant at early times, in other words 
that the angular size of the peaks of the CMB and baryon oscillations would be the 
same as in ACDM. Later we will go on to explicitly show that this is a bad assumption. 

The potential U(a) for the f(R) = R — /// R model is plotted in the first diagram 
of figure HI it rises from zero to a maximum and then falls back to zero as a goes to 
infinity. 

Figure [1] shows the combined best fit regions for the supernova, BAO and CMB 
data for Qm = 0.3 plotted as a function of o"o and da/dr\ Q . The best fit to supernova 
data only is obtained for o"o = 0.22 and da/dr\ = —0.24, and the best fit when also the 
two angular distance data points are included for a = 3.20 and da/dr\ = 1.16. Plotted 
also are look back times for the age of the universe in this model assuming a Hubble 
constant H = 65kms~ 1 Mpc _1 (or more precisely the elapsed time in the matter frame 
since a redshift of z — 20) and a line labeled "big bang" which divides the parameter 
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Figure 1. Look back times (age of the universe) for f(R) = R — fJr/R cosmologies. 
We label different cosmologies using the value of the effective scalar a and its time 
derivative in the Einstein frame da/dr at redshift z = 0. These two parameters plus 
the matter density (this plot is for £Im = 0.3) provide three of the four parameters 
necessary to determine the cosmology. The fourth is the Hubble constant which is 
integrated over in cosmological fits. Look-back times are in years in the matter frame 
assuming a Hubble constant of 65 kms _1 Mpc _1 . In cosmologies above the line labeled 
'big bang', the scalar field goes over the top of the potential too quickly (see text). 
The red region is the 99% region fitting to supernova data, the blue region is the same 
for BAO. 

space between the grey region of solutions where the scalar field a — > too quickly in 
the past, which corresponds to a curvature singularity 

This can be understood in the following way, for a given value of cq, increasing 
values of da/dr push a higher up the effective potential as one moves into the past, 
which means that the rapidly redshifting contribution from Qm can be reduced and 
the universe expands more slowly in the past, making it older. This increase in age 
with do~/dr\ ends abruptly for some value of da/dr\ where a goes over the top of the 
potential and falls to o = relatively recently in cosmological terms, which signals that 
R — > oo in the matter frame. That region beyond the curve labeled with the words 
"big-bang" is therefore not included in the analysis and x 2 values are not calculated 
there. 

Also plotted are the 99% confidence bands for the supernova data in red and the 
baryon oscillation data in blue. These regions overlap each other for certain values of 
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Qm and therefore seem compatible with the model. 

The best fit \ 2 values are all rather close to the dividing line to the region where 
these maximum age cosmologies exist. In these models, the scalar field runs back from 
the top of the potential so that a — > oo, R — > and the effective total equation of state 
w — > 1/3. Speaking purely in terms of expansion history, the universe therefore becomes 
effectively radiation dominated very recently as we go back in time, not because of the 
energy density of any radiation, but simply because of the modified gravity giving rise 
to what looks like, from the FRW perspective, an energy density with equation of state 
w = 1/3. This has been noticed by other authors [301 13"T] . 

On the right of figure [2] is plotted the effective total equation of state that one would 
obtain if the same expansion was due to some energy density rather than a different 
theory of gravity, in other words 

2 (1 + *) dH 

W= 3^-^~ 1 (16) 
and it can be seen that at high redshifts, w — > 1/3 which is the same as pure radiation. 

On the left of figure [2] we can see that H(z)/H is similar in both cases at low 
redshifts where the supernova and baryon oscillation data is fitted, then the H(z) for 
the f(R) model dips below, then rises above, the H(z) for ACDM. This means that the 
luminosity/angular distance integrals between z = and z = 1000 are similar. This 
point is made clearer by looking at figure [3] which shows the fractional difference between 
the luminosity distancj§| for the f(R) model and the ACDM model. The difference 
between the two is rather small at very low and high redshifts but peaks around z ~ 5. 

Better knowledge of the Hubble diagram around redshifts of z ~ 1 — 5 would 
therefore probably differentiate between the models, but since we only have good data at 



§ and therefore also the angular distance since they are related by a factor of 1 + z 
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Figure 3. Difference between the luminosity distance as a function of redshift for fiat 
ACDM and the best fit 1/R model when VLm is taken to be 0.3. 

z < 1.5 from supernovae and at z ~ 1100 from the CMB we are not able to differentiate 
between them using only these data sets. 

Figure H] shows the comparison of the data with the f(R) = R — fi 4 /R model for 
different values of Qm, &o and da/dr\ . Plotted are the banded constraints corresponding 
to the supernova data and to the baryon oscillation data. A large region of the 
parameter space is ruled out as it corresponds to regions where o —>■ too quickly 
in the past, corresponding to a curvature singularity which would be too recent to 
accommodate the early universe physics that we know must take place (last scattering 
surface, nucleosynthesis etc.). The best fit values of x 2 f° r different values of Qm are 
listed in tableland they show that the model works well in obtaining expansion histories 
rather similar to ACDM. Addition of the CMB data makes the fit worse but it is still 
possible to fit the data with a \ 2 P er degree of freedom which is less than one. Also the 
look back time is consistent with the age of globular clusters [32] . 

The same analysis has been carried out for the f(R) = R — /U 6 / R 2 model of gravity 
and the results are listed in figure [5] and the best fit values of x 2 are listed in table [U 
The situation is completely analogous to the —fi 4 /R modifications. Again, the models 
are able to explain the supernova data. 

To summarize this section, universes which fit the data in the 1/R model lead to 
the scale factor a oc t 1 ! 2 at high redshift. Despite this, the universes obtained are old 
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a. 



Figure 4. Comparison of f(R) = R — fi^/R with data. The top left hand panel is the 
effective potential in the Einstein frame, while the other three plots are the values of 
(Jo and da/dr |o favoured by the supernova data in red and the BAO data in blue (see 
caption of figure [T] for more details explaining the plot). The three plots correspond to 
different matter densities - CIm = 0.2, 0.25, 0.3. 
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xLn(SNe+BAO+flat) 


169 


156 
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156 


Xmin(SNe+BAO+CMB+flat) 


201 


182 


179 


186 


164 


159 



Table 1. The best fit values of x 2 (178 d.o.f. when the supernova data only is included) 
for the —fjL 4 /R and the —fi 6 /R 2 modifications of the Einstein-Hilbert Lagrangian. The 
vertical columns correspond to different values of J7m while each row corresponds to a 
different subset of the data. Flatness is assumed throughout. 
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Figure 5. Comparison of f(R) = R — /J. 6 /R 2 with data. The top left hand panel is the 
effective potential in the Einstein frame, while the other three plots are the values of 
<7o and da/dr |o favoured by the supernova data in red and the BAO data in blue (see 
caption of figure [T] for more details explaining the plot). The three plots correspond to 
different matter densities - Qm = 0.2, 0.25, 0.3. 

enough to accommodates globular clusters. Also, because the integrals over 1/H(z) are 
dominated at low redshifts, it is possible to obtain expansion histories at low redshifts 
rather similar to ACDM. 

As we mentioned at the beginning of this section, this ability to fit data which is 
relevant for the ACDM universe is not at all interesting unless the data itself is still 
relevant. While this is true for the supernova data, the different behaviour of the scale 
factor at early times in these inverse R models will change the CMB, and therefore 
trying to fit the expected position of the first peak that one would expect if the early 
universe contained only matter and radiation as we have done above would be a mistake. 
The simplest way to see how the different expansion in the early universe would affect 
the CMB is by looking at the thickness of the last scattering surface. 
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4. Thickness of the Last Scattering Surface. 

So far it has been shown that in models of modified gravity of the form 1/R and 1/R 2 , 
the expansion of the universe can fit the supernova data. 

The age of the universe in these models is perfectly compatible with the age of 
globular clusters observed. One therefore cannot use such naive benchmark tests to rule 
out the model. The integral of the inverse of H(z) between z = and z = 1000 in these 
models is not too different from ACDM so that features with the same physical size and 
redshifts as the peaks in the CMB and the galaxy correlation function would posses the 
same angular size in this universe. However, this is not encouraging since the physical 
size baryon oscillations and also the peak in the CMB would be completely different due 
to different early universe physics. 

There are of course still unused weapons in our arsenal of observations, in 
particular one can think about how the perturbations grow in these models, as has 
been analysed previously in references [33l IMl [19]. Also one can consider what happens 
to nucleosynthesis which is in principle rather sensitive to the expansion of the universe 
during the freeze out of the weak interactions. 

The size of the horizon at the last scattering surface will also be different due to 
the different expansion which occured at higher redshifts and consequently the time 
corresponding to that epoch. This would directly shift the angular position of the 
acoustic peaks. Furthermore structure formation would be different, leading to very 
different predictions for the size of the peaks. 

However, it is simpler and easier to consider the epoch of recombination and to 
look at the way that the added expansion at redshifts of order z ~ 1000 would affect 
the thickness of the CMB. 

The re-combination of protons and electrons which occurs as the universe cools 
is the process responsible for making the universe transparent. In normal cosmologies 
where the universe is still matter dominated at this epoch, the electrons, photons and 
protons remain rather close to thermal equilibrium throughout the process. The physics 
is therefore less dependent upon the expansion rate of the universe but more upon the 
temperature. However, the effective radiation density that one obtains at high redshifts 
for the f(R) models considered here is orders of magnitude larger than the energy 
density of radiation implied from the temperature of the CMB. 

The density of radiation implied by looking at the 2.7 K background radiation 
is of the order of fi 7 ~ 10~ 4 so that matter radiation equality occurs at redshifts 
around a few times 10 4 . The effective radiation density that fits the H(z) at high 
redshift correspondent to the best fit 1/R model (to the supernova data alone) is 
n„ =1 / 3 ~ 0.22. The expansion rate around z ~ 1000 is therefore changed considerably, 
and recombination takes longer to occur. 

The equation which needs to be solved in order to calculate the rate of reionisation 
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is j33 Eg 



dx P 1 



dz H(l + z) 



an p x e — (3(1 — x e ) exp 



B\ — B-2 
kT 



C (17) 



where a is the recombination coefficient, (3 is the ionization coefficient and B n is the 
binding energy of the n— th level of the hydrogen atom. The factor C is given by 

l + KT(l-x e ) 

l + K(T + (3)(l-x e ) { } 

where T is the decay rate of the 2s excited state to the ground state via the emission of 
2 photons. The ionisation coefficient (3 is given by 

f 27rm e kT \ 3/2 f B 2 \ 

and the recombination coefficient a is described by the expression 

a _ E t (2f + l)87r / kT V /2 (B n \ r QniV 2 dy 
c 2 \2nmJ \kTJ J Bn/kT exp(y) - 1 

Solution of equation (1171) gives the ionisation fraction as a function of redshift. Following 
|35j . we approximate a to be a oc T~ a5 and normalise it so that it gives rise to the 
correct redshift for ionization as observed by WMAP. The thickness of the CMB Last 
Scattering Surface (LSS) can then be found by looking at the probability for finding a 
photon reaching us from a redshift z, g(z) = e~ T dr/dz. We then define the thickness 
as being the full width at half the maximum of this function g(z). In this way, when 
we consider recombination in the ACDM universe we obtain a thickness for the CMB 
of Az = 197 rather close to the WMAP value Az = 195. 

Figure [6] shows us the thickness of the CMB as a function of the effective radiation 
density. The greater the expansion of the universe at last scattering, the more gradual 
is the era of recombination (gradual with respect to redshift rather than time) and the 
larger the thickness of the LSS. The redshift of the LSS also changes, but the fractional 
change in the thickness is much greater. The thickness derived from observations by 
WMAP is Az = 195 ± 2 [37], however we find that for the best fit parameters to the 
supernova data alone, we get a thickness much closer to Az = 530. 

The first effect of a shift of the last scattering surface to earlier times would be a 
different sound horizon size at decoupling, so that the peaks of the CMB would change 
position [38]. At the same time, the increase in the thickness of the LSS would increase 
photon diffusion within the surface, blurring out features with a physical size much 
smaller than the LSS. Suppression at high multipoles would therefore be brought down 
to much lower multipoles [39] . 

Constraints upon the total number of relativistic species present which are not 
chemically coupled to the plasma have been obtained and refined by running codes like 
CMBFAST which take into account both the effect of the shift in the peak and the 
supression of higher multipoles (see [3H] and references therein). Those studies show 
that an increase in the number of relativistic degrees of freedom of factor AN = 13.37 
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Figure 6. Thickness of the LSS as a function of the effective radiation like energy 
density. The red line runs from the value predicted from the CMB in ACDM to the 
value which best fits the high redshift expansion for the 1 /R model which best fits the 
data. 



can be ruled out at 2-cr by looking at the resulting fit to the CMB data alone. Since 
the modified gravity solutions which fit the data found in this paper correspond to an 
increase in the effective relativistic degrees of around AN ~ 300 at high redshift, one 
can clearly rule them out. 

In summary, although the 1/R models fit the existing supernova data, they can be 
ruled out because they will lead to a LSS which is too thick. The theory therefore is 
inconsistent with the data it is designed to fit. 

4-1. Finding a solution which fits with the CMB 

We have seen that the solutions which fit the supernova data best lead to problems in 
the early universe due to the extra radiation-like effective energy density which occurs 
at that time. An important question is whether or not it is possible to find solutions 
to the equations where the expansion rate at the redshift corresponding to the last 
scattering surface is comparable to the one predicted in ACDM cosmology, which is in 
good agreement with the observed rate of reionisation. 

We have seen that the thickness of the LSS is increased to Az ~ 530 for the best 
fit to the supernova data alone. If we try to use the CMB data and the BAO data 
as angular sizes to fit also, the thickness of the LSS in the best fit case decreases to 
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Figure 7. Comparison of the best fit version of f(R) = R — fJ. A /R (red curve) and 
ACDM (green curve) vs. a f(R) — R — H 4 /R solution that is close to the ACDM at 
the redshift of decoupling (black curve), in all cases Qm = 0.3. 

Az ~ 350. This is still far to large to make the use of the CMB and BAO angular sizes 
valid, but it shows that we need to check to see if it is possible to obtain parameters 
where the thickness of the LSS might be consistent with the CMB data. 

Figure [7J shows the H(z)/H for ACDM, for the best fit 1/R model mentioned in 
section [7] and a third cosmology. This third cosmology is a 1/R model where the Hubble 
rate at reionisation is close to that experienced in ACDM as can be seen in the figure. 

The x 2 fit to the supernova data for these parameters is 184.5, which, while much 
worse than the best fit model described earlier, or the best fit to ACDM, is still close to 
unity per degree of freedom. 

What also can be seen is that at lower redshifts, H(z) is much lower than 
the corresponding expansion rate in the ACDM universe. In fact the age of the 
universe where the expansion rate at the last scattering surface fits with observations 
is t — l-6/i?o which corresponds to 25 Gyrs for h =0.7. In such a universe one would 
expect to find globular clusters with a much lower turn off in the HR diagram then what 
is observed, the current inferred result for the age of the oldest clusters being 13.2 ±1.5 
Gyr jH]. The white dwarf cooling sequence in globular clusters also points to similar 
age estimates, making it difficult to live with such and old universe [42J. 
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In this paper we have put cosmological constraints on models where the dark energy 
component of the universe is explained by theories where gravity is modified by the 
addition of terms to the Einstein Hilbert Lagrangian which contain inverse powers of 
the Ricci scalar R. 

We find that both 1/R and 1/R 2 models, when forced to fit the supernova data, 
give rise to solutions at high redshifts where the scale factor a oc t 1//2 . The universe 
does not therefore have a matter dominated phase, but rather interpolates between an 
accelerating phase at low redshifts and a 'radiation-dominated' phase at high redshifts. 
The expansion at this earlier phase is driven not by radiation which we know the density 
of by observing the temperature of the CMB. Rather this radiation-like expansion is due 
to the modified Friedman equations in this class of models. (More complicated models 
may be able to avoid this feature, see |40j.) 

Despite this, we can still find solutions for the Hubble expansion H(z) that are 
consistent with the recent supernova results. We have however shown that the theory 
leads to very different predictions as to what the CMB and therefore the baryon 
oscillation data should be if it is to fit the luminosity and angular distance data between 
redshift zero and the last scattering surface. If we try and force the theory to fit the 
supernova data while obtaining the correct thickness for the last scattering surface, we 
are then left with a universe which is much older than the oldest stars observed in our 
own universe. We therefore argue that these models are inconsistent with the data. 
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